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Making use of inverse Mellin transform techniques for analytical continuation, an 
elegant proof and an extension of the zeta function regularization theorem is obtained. 
No series commutations are involved in the procedure; nevertheless the result is naturally 
split into the same three contributions of very different nature, i.e. the series of Riemann 

43 ■ 

zeta functions and the power and negative exponentially behaved functions, respectively, 
well known from the original proof. The new theorem deals equally well with elliptic 
differential operators whose spectrum is not explicitly known. Rigorous results on the 
asymptoticity of the outcoming series are given, together with some specific examples. 
Exact analytical formulas, simplifying approximations and numerical estimates for the last 
of the three contributions (the most difficult to handle) are obtained. As an application 
of the method, the summation of the series which appear in the analytic computation (for 
different ranges of temperature) of the partition function of the string — basic in order to 
ascertain if QCD is some limit of a string theory — is performed. 
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1 Introduction 



Zeta-function resummation formulas are essential in the zeta-function regularization pro- 
cedure . In particular, they constitute the key ingredient for the proof (done in several 
steps by Weldon, Actor, Elizalde and Romeo 0, ||) of the zeta function regularization 
theorem. The importance of the final outcome of the theorem — and also, of its extension 
to multiple, generalized Epstein-Hurwitz series with arbitrary exponents || — has been 
stressed by several authors Q. Specifically, it is necessary for actual, practical applications 
to estimate or put a bound to the error one introduces by neglecting the non-polynomical 
terms which arise in the series commutation process. These terms have proven to be 
very difficult to handle and, due to the fact that they are (very) small in some cases, 
usually they have just been dropped off from the final formulas. To put an end to such 
an unpleasant situation is one of the purposes of this paper. 

Another motivation is to present a completely different method for the derivation of 
these additional contributions, which comes about from a very elegant approach that has 
its roots in an elaborate admixture of the Mellin transform technique and the heat kernel 
method. Our new procedure yields a very convenient, closed expression for the zeta 
function corresponding to very general elliptic operators in terms of complex integrals 
over movable vertical lines in the complex plane. The final result is the same as the 
one painstakingly obtained via the original method of the series commutation — with the 
great advantage that the proof of asymptoticity of the resulting series turns to be now 
immediate and also, that the calculation of bounds to the additional term is considerably 
easied. Moreover, cast in this form the theorem admits a quite natural generalization to 
the case where the spectrum of the operator is not known. 

In the next section of the paper, we describe the Mellin transform method as ap- 
plied to generalized heat kernel operators. The method is extended in section 3 to the 
case of trace formulas of generalized zeta-function type. In section 4 some mathematical 
theorems on asymptoticity of series, basic for practical applications of the procedure in 
rather general situations are presented. In section 5 we obtain some initial (but already 
non-trivial) results on the additional term by making convenient use of the Poisson sum- 
mation formula. Specific numerical evaluations and an analytical bound on the elusive 
additional term which appears as a byproduct of, e.g., the series commutation procedure, 
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are obtained in section 6, together with some convenient integral expressions for this term. 
As a particular example of the usefulness of the whole procedure, in section 7 we obtain 
exact expressions for the analytical continuation of highly non-trivial, zeta-type series 
which appear in the calculation of the partition functions of strings and membranes. This 
issue happens to be essential in deciding if QCD can (or cannot) in fact be described as 
a (super) string (or membrane) theory. Finally, section 8 is devoted to conclusions. 

2 Inverse Mellin transform method for generalized 
heat kernel operators 

We shall here rederive the (by now well known) zeta-function resummation formulas 
— using the Mellin transform techniques that have been employed for obtaining high- 
temperature expansions, both in field theory as well as in string theory || — in a very 
straightforward way, 

Let us consider a non-negative elliptic operator A of order d, defined on a compact n- 
dimensional manifold. Let {A«, <pi\ be the spectral resolution of A. (The explicit knowledge 
of the spectral resolution of the operator will not be basic to our method, and we shall 
make only a very limited, circumstancial use of it.) In what follows we shall be interested 
in the quantity 



which can be interpreted as a generalized heat kernel operator. It satifies the functional 
differential equation 






(2) 



and, from the operator identity 




(3) 



one obtains the complex integral representation 




(4) 
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Since A is assumed to be elliptic, the diagonal part of its heat kernel as well as the 
kernel of the related zeta function exist. A theorem of Seeley yields 



( A (s + z)(x) 



T(s + z) 



E 



KJx) 



S + Z + 



+ J(s + z, x) 



(5) 



(where ~ is the symbol for asymptotical equivalence), being K r (x) the heat kernel expan- 
sion coefficients and J the analytical part of (a, which is not known explicitly, in general. 
As a result of all the preceding expressions, one arrives at 



F t (s,x) ~ 



c T(s + z) 



E 



KJx) 



S + Z + 



+ J(s + z,x) 



2ni JRez 

The corresponding trace formula can be obtained by integrating over x, with the result 

K 



(6) 



Tr F t (s) ~ Yl K Se ~ tX ' = /j 



Rez = c l(s + z) 



E 



r s + z + d 



- + J(s + z) 



(7) 

In the rest of the section we will employ this form of complex integral representation, 
with the understanding that all the conclusion we shall draw using this method can be 
also validated for the corresponding non-integrated quantities. In order to obtain the so- 
called resummation formulas (see Actor, Weldon, Elizalde and Romeo [fj, |3|) one needs to 
perfom a very accurate application of the residues theorem. For a generic s, we encounter 
simple poles at z = —k, k = 0,1,2,..., and at z = " n -^- — s in the expression above. 
Furthermore, in general there will be a non-zero contribution coming from the integral 
evaluated along the semicircumference at infinity (see Elizalde and Romeo, refs. [§, [H). 
As a result, we can decompose F t (s) into three parts, very different in nature, 



Ft(a) 



~Fi 1) (s)+Fr(s)+Fyis ) 



,(2) ( 



,(3), 



which turn out to be, in our language, 



F t (1 \s) = j: { -^rUs-k), 

k 



F t (2 \s) 



and 



(9) 
(10) 

(11) 



C being a convenient contour on the complex plane. Some remarks are in order. First, 
these results are in complete agreement with the ones originally obtained by one of the 
authors (see also the recent paper by Actor ||) for specific situations. However, in our 
new approach, we explicitly see how the heat kernel coefficients K r enter the game, in the 
general case considered from the begining (eq. ([!])). These are, in principle, computable 
quantities (for a recent reference see Branson and Gilkey 0), in terms of which the 
characteristics of an elliptic operator are always given. The analytic part of the final 
result — which is generally unknown — enters into the first and third contributions. In 
some cases (see, for instance, Elizalde and Romeo M), it can be explicitly evaluated. In 
particular, it is interesting to consider also the case of hyperbolic compact manifolds, 
where the analytic part can be related to a zeta function of Selberg type (Bytsenko et 
al. H), even though yet in this case the spectrum is usually not known. Also, we should 
notice that a similar result has been obtained by Avramidi in a particular case ||. 

A very important remark concerns the range of the sums in eqs. (9)-(ll). In previous 
works |2], |^, |J this range was formally taken to be infinite: the sums were series and 
the integration contour C in was at infinite distance from the origin. Then these 
expressions can only have a formal character: on the one hand the series converge only 
in very specific cases and, on the other, the contribution of the contour is usually infinite. 
This will not be any more our philosophy in the present paper. In the following sections 
we shall prove that (under very general conditions) the series in (9)-(ll) are asymptotic; 
thus, we will always consider them as cut after the most favourable term (yielding always 
finite sums), and the contour C will be a vertical line at a constant finite abscissa c (i.e., 
Re z = c), together with the two corresponding horizontal segments at infinity. This is 
the central part of the full method in this paper, and must be kept always in mind, from 
now on, when considering expressions like (9)-(ll). Rigorous specifications on this point 
will be given in the next sections. 
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3 Extension to trace formulas of generalized zeta- 
function type 



The above technique can be extended to trace formulas of the kind 

F(s,t) = J2K s f(t\i), (12) 

i 

provided f(y) admits an invertible Mellin transform, namely 

f(y) = ^-f dzy- z M[f](z), (13) 

2m J Re z = c 

with 

M[f](z) = / dyy^fiy). (14) 
Jo 

Being more precise, it is assumed that the function / is integrable over the open 
interval (ti,t 2 ) provided that < t\ < t 2 . Defining 

a = M{a* eR\f{0 + ) = O{r a *)}, j3 = sup G R | < /*(+oo)0(t~ /J *)} , (15) 

we shall also consider that a < (3. This is no restriction since if, on the contrary, — oo < 
P < a < oo, one just needs to define (for R > 0) j\{t) = 6{R-t)f{t), f 2 {t) = 9{t-R)f{t) 
and Fj(s, t) = J2i ^7 s fj(t^i), j = 1, 2, and then consider each Fj separately (with ai = a, 
fti = +oo, and a 2 = — oo, (3 2 = +/?)• 

In general, there will exist a value s G R such that for Re s > s the series J2i 
will be absolutely convergent and for Re s < so it will be not. Then, a sufficient condition 
in order that ( |T2| ) makes sense is Re s + (3 > sq. In this case, 

with c G (a, (3) and Re (s) + c > so- We take M[f](c + iy) G L 1 (— oo < y < oo), so that, 
finally, making use of all these equations, we can write 

F(s,t) = ^- [ dzt- z M[f}(z)C A (s + z). (17) 

2m JRez = c 

Now, we must suppose that one can solve the problem of the analytic continuation of 
the functions (a( z ) an d M[f](z) to the left in Re z. Under very general conditions, the 
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function M[f] can be extended as a meromorphic function, but for a general (a this is not 
so easy to do in practice. In general one must calculate all the Seeley-De Witt coefficients 
of A and use them to infer the meromorphic behavior of (a (in the second example below 
we show an alternative way to extend (a using our knowledge of the Aj only). We also 
need to know the behavior of (a(z) and M[f](z) for |Im z\ — > oo. 
As before, for a generic s, we can split 



where now 



F(s,t) = FW(s,t) + F®(s,t) + F®(s,t), 



fM(s, t) = ResM[/](z fc ) t- z % A {s + z k ) 

r 1 { d > 



(19) 
(20) 



and 



F^(s,t) = ^- I dzM[f}(z)r z Us + z), (21) 

Z7TZ JC 

being z k all the (simple) poles of M[f](z). If double poles occur, logarithmic terms in t 
show up. This may happen for particular values of s, once the other free parameters are 
held fixed. 

Let us now consider some illustrative examples. 
Example 1. For the Riemann zeta function itself, we have 



C(s) = C(|t| p(r) ln|t] 



s = r + it, 



(22) 



with 



p(r) 



1 

2 


r, 


—K <r< -6, 


1 
2 


r, 


-5 < r < 0, 


1 
2 


r, 


< r < 5, 


1 - 


r, 


5 < r < 1 - 5, 


1, 




1 - S < r < 1, 


0. 




1 < r < 1 + 5, 


0. 




r>l + S, 



(23) 



where K > and 5 > (a standard choice is S < 1 and K arbitrarily big). Actually, In \ t\ 
can be supressed in the first, fourth and seventh intervals. It is remarkable that eq. ([22]) 
is valid uniformily in each of the intervals of (f23|) . 
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Example 2. Take A n = n and add a degeneration g(n) (being g an analytical function 
with good behavior in order to fulfill what follows), i.e. 

<.W-f ^ = ~/ c M^)]|(r*,)*, P4) 

being C the contour of the sector of the complex plane defined by | argz| < # and \z\ > e, 
< e < 1. Decomposing C into its upper and lower parts, C*, and assuming g to be 
analytical on the sector, we obtain 

roc A j if/ \ d 

In general, the second and third terms on the r.h.s. are integer functions of z, and only 
the first one needs to be continued to the left in Re z. This can be done by applying the 
Mellin transform techniques. 

4 A more explicit study of the additional term 

Of the three terms which appear in each of the splittings (S) and (|T8|), the third one is 
the most difficult to tackle. In fact, it was completely overlooked in the first formulations 
of the zeta function regularization theorem, what led to several erroneous results. In the 
procedure above it has shown up in a nice and elegant way. However, its apparently very 
simple form is somewhat deceiving. In fact, in its primitive form P, [5], it required an 
evaluation of the integrand — which involves the zeta function of the elliptic operator — 
all the way along a contour of infinite radius. We shall here show explicitly how to obtain 
an alternative asymptotic expansion of expressions (H) and (|T8|) in which the residual term 
consists of a complex integration on a vertical line of constant, finite abscissa. We start 
by proving the following theorem. 

Theorem 1. Let 6 > and h be an analytic function in the domain defined by z ^ 
and | argz| < 9 , and such that 

oo 

h(z)~5> n A |*|->0, (26) 

n=0 

and (only to simplify the discussion, by putting /3 = oo) 

h(z) < O (z~ R ) , V.R G R , \z\ -»• oo, (27) 
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where we understand that Re (3 n < Re (3 n+ i and Re (3 n — > oo. Define 

h{n r r) 



n=l 



rr 



(2* 



with r, r > and s G C. Then, Y r (s,r) is an integer function with respect to the 
argument s and, denoting as before the Mellin transform of a function / by M[/], one 
has 



YJs,t) ~ T 



-(l-s)/r 



-M[h) f^-i 
r V r 



M[h]' ( — 4 I 



1 



.s/r 



'-M[h] 



M[h]' 1 



(29) 



n=0 



Proof. Write 



n(s,r) 



1 r _ 1 r d r _ i 

/ dz 7r cot(nz) z~ s h(z r r) = / dz In sm(irz) — z~ s h(z r T) , (30) 

2ni Jc 2ni Jc dz L i 



C being the circuit made up of the two straight lines at angles ±#o- From the first 
hypothesis it turns out that the two integrals in which ([$(]) splits exist for any z ^ with 
I argz| < 6 — 5, V5 > 0. Thus we deform the circuit C to C\ corresponding to the lines 
at angles ±(#o — $) an d with a small circular arc of radius e at the origin. 
Notice now the asymptotic behavior of a function of the following type 



T-Ci-s)/ 1 " /I — S 

dpp- s h(p r T) M[h] 



S a ™ r g , i=£ ' 

n=0 ' «T r 



(31) 



provided there are no pole coincidences. Then, defining 

Y 1 = — ( dz lnsin(Trz) z-'^hifr) = Y+ + Y{ 
2lT Jc 



(32) 



corresponding naturally to the two branches of the circuit, C±, the analysis of each of 
these integrals is quite simple. The second part can in fact be treated as the first, by 
noticing that 



r r 



Y 2 ee / dz lnsin(Trz) z r -'- L h'(z r T) 

2iri Jc 



—YSr h \s-r,T) 



s — r 



(33) 



Collecting everything together, one arrives at the desired expression fl29|) , valid in prin- 
ciple for any positive integer n such that s ^ 1 + /5„r and s ^ (3 n r. These cases must 
be considered specially. Either one starts over again or, what is better, one proceeds 
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according to the following hint: to perform a shift s = 1 + f3 n r + 5 or, correspondingly, 
s = f3 n r + 5, and then make use of the well known lemma 

Lemma 1. Let D an open domain of C, s G D, D* = D — {s }, and f3 n a sequence 
in C such that Re /3 n < Re f3 n+ \ and Re f3 n — ^ oo. Let Y(s, r) analytic as a function of 
the first argument for s G D and r G (0, cc), a > 0. Suppose that for Vs G -D* one has 
Y(s,t) ~ X^o a ™( s ) r/3n 5 where the a n (s) are analytic on D, and that this expansion is 
valid uniformily for s G C, C being a path contained in D and enclosing s . 

Then it turns out that also Y(s , r) ~ J2^=o a n (s ) T^ n . The proof is very simple. 

In order to determine the behavior of M[f](x + iy) for \y\ — > oo the following theorems 
are useful 

Theorem 2. Let / G C n (0, oo) and supose there exists x G R such that, for re > x , 
lim t _ >0O (id/ 'dt) p (t x f) = 0, for p = 0, l,2,...,n, and that t _1 / p is absolutely integrable, 
where f p (t,x) = {td / dtf {t x f) . Then M[f](z) = 0{\y\- n ) for y -> ±oo. 
Proof. Follows immediately by wishful application of Riemann-Lebesgue's lemma. 

Theorem 3. Let / G C n (0, oo) with the asymptotic behavior for t — > + : /(i) ~ 
Y^m=oPmt am , with Re a m monotonically increasing towards infinity. Suppose also that 
asymptotic expansions for t — > + for the succesive derivatives of / are obtained by 
taking the derivative, term by term, of the above expansion for /. Suppose also that 
lim^oo f p (t, x) = 0, for p = 0,1,2, ... ,n, and that t~ x j v is absolutely integrable for 
x > — Re a . Then it follows that M[f](z) = 0(\y\~ n ) for y — > ±oo and for any x, that 
is, in the whole domain of the function, after its continuation. 

Proof. Let p G R be large enough and /i such that Re a M _i < p < Re a M , where for a pos- 
itive integer 5 satisfying Re a + 5 > Re a M , we define <J p (t) = exp(— t 5 ) ^m=oPmi' lm , and 
/ = f—Op. Applying the Theorem 2 to / and noticing that M[a] p (z) = Y^!^} (p m /5)r((z+ 
a m )/8) has exponential behavior for \y\ — > oo, from M[f] = Mf + M[a] p (which gives the 
continuation to -Re a p < Re z < (3) and the fact that p is arbitrarily big we prove the 
theorem. 

For xq, 6 G R, let us define that a function h belongs to K, h G K(xq, 9q) iff for any 
x > xq and e > we have M[h](x + iy) = O[exp(—(9 — e)\y\)], for \y\ — > oo. Also, denote 
the sector S(9 ) = {t ^ 0| | argt| < 9 }. 
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Theorem 4. Supose that in the sector S(9 ): (i) h is analytic; (ii) h = 0(t a ), t — > + ; 
and (iii) h = exp{-at u ) E~ = oES ) c mn (lnt) n )T r ™, t -> oo, with Re a > 0, v > and Re 
r m monotonically increasing towards infinity. Then: (i) if a — 0, h & K(— lie a, do); (ii) if 
Re a > 0, h 6 Re a, 0), where 6 1 = min(6 l , (-7T - 2| arga|)/(2i/)). 
Proof. We only do the first case (the second is analogous). Observe that M[h](z) is 
analytical for Re z > — Re a. Define 9' = 9 — e, 9 as given in the theorem and e > 0; 
deforming the circuit by an angle ±9' and making the change of variable t = re ±l6 ' (r, 
real, will be the new variable of integration) in the Mellin transform integral, one easily 
shows that M[h](z) = 0[exj>(—9'\y\)), for \y\ — > oo, x > Re —a. 

With a similar strategy as that used in Theorem 3, with convenient choices of the 
function a p , one can prove results analogous to Theorem 4 (so called theorems of expo- 
nential decrease) valid for a range of x which moves towards the left. The only proviso is 
that h have a convenient asymptotic expansion around 0. 

Also to be noticed is the fact that in order to apply these theorems, with the aim of 
displacing the contour of integration in ([T7D, one must make sure that the conclusions 
of Theorems 2-4 are valid uniformily on any segment of the x variable along which we 
want to displace the integration contour. If we can perform the translation from x = c to 
x = c', d < c, and we check that (a(s + d + iy)M[f](d + iy) is absolutely integrable, then 

F(s,t) = J2 Res [Us + z)t' z M[f}{z)\ + l — J Rez = d dzUs + z)t'^M[f]{z). 
c'<Rez<c 

(34) 

To finish this section, let us consider a different example which cannot be resolved by 
direct application of Theorem 4 above (since it corresponds to a ^ and Re a = 0). 

Example 3. Let 

F( M ) = £^M, (35) 

n=l " 

that is, X n — n, f — Jfi, a — —fi, (3 = 1/2. We have 



2 



with So = 1, Re s > 1/2 (in this case the domain of s can be enlarged), and c > 1— Re 
s < 1/2. M[J}^ has poles for z = —fi+2n, n 6 N , with residues (-l) n 2~^~ 2n /[n! (/i+n)!], 
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while ((s + z) has a pole at z = — s + 1. From (|3T)| ) and from the behavior of the gamma 
function we know that in any closed interval of x = Re z we get the behavior 

M[Jl(z) = 0(\yr^), Ij/I-oo, (37) 

uniformily in x. We have Re s = 1/2 + ex, ei > 0, and choose c = 1/2 — 62, with 
< e 2 < ei, and any 5, < 5 < min{l/2, E\ — e 2 }. 

Consider now the following, segmentwise, displacement 

G(z)=M[JUz)((s + z). (38) 

Then it turns out that 

1. for 1/2 + 5 - e x < x < 1/2 - e 2 , \M[J]^z)\ = 0{\y\- 1 ^) and |C(s + z)\ = 0(1), 
therefore 

|G(z)| = 0(|yr 1 / a -«; (39) 

2. for 1/2 -e a < x < l/2 + «J-ei, |M[J] M (*)| = O^' 1 ' 2 ^) and |C(s + *)l = 0{\n\y\), 
therefore 

IG^^Odyl-^lnlyl); (40) 

3. for -1/2 + S - ex < x <l/2-ex, \M[J]^{z)\ = 0{\y\-^ 2+& -^) and \({s + z)| = 
0(|2/| 1-i ln|2/|), therefore 

IG^^Odyl-^lnlyl); (41) 

4. for -1/2 - e x < x < -1/2 + 8 - e u |M[J] M (z)| = 0(\y\-V 2+5 - ei ) and |C(s + z)\ = 
Odyl 1 / 2 ), therefore 

\G(z)\ = 0(\y\-^); (42) 

5. finally for -K-l/2-ex <x< -1/2-ex, \M[J}^z)\ = O^y^' 1 ) and |C(s + *)| = 
(D(\y\~ x ~ ei ), therefore 

\G(z)\ = 0(\y\- 1 -«), (43) 

where K is arbitrarily large. 
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5 The additional term and the Poisson resummation 
formula 



We would like first to examine a particular but important case, as an illustration of the 
difficulties associated with the determination of the additional term. Such case will be 
considered also in the next section but under a different point of view. 

Let us consider the selfadjoint operator A = H@, where H is the Dirichlet Laplacian 
on [0, 1] with eigenvalues A n = n 2 . Putting 2(3 = a, we can write 



( A (z) = Y,n- az = aaz). (44) 

n=l 

The related heat kernel trace reads (t > 0) 

Tre" M = = S a (t). (45) 



1 r r(-) 00 (— t) k 

Sa(t) = — / dzt-*C(az)T(z) = + Y / - 1 J-a-uk) + A a (t) 

ZTXl JC fc! 



Making use of the Mellin representation discussed in Sec. 2, we get 

at a k=0 

T(-) 1 00 (-t) k 

ata A k=l K - 

where the power series in t will converge for \t\ < b (for some b, abscissa of convergence). 
Then we can analytically continue to the remaining values of t. The quantity S a (t) may 
be evaluated by making use of the Poisson formula, which states that for f(x) = f(—x) 
and f(x) G L\, the following equation holds: 



00 j roo 00 roc 

J2f( n ) = -ofW+ dxf(x) + 2Y J dxf(x) C os2nnx. 



(47) 



Let us consider the function 

/(x) = e-*N a . (4 
An elementary computation permits to write 

at 



1 T(-) 00 roc 

S a (t) = -- + -^f- + 2 E / dxe-^ cos2imx 

2 n f~ n=l Jo 



(49) 



and by comparing eqs. ( |46f) and fl4"9"|) , we get 

{-t)\ 



^ poo ^ + * 

A a (t)=2j2 dxe- lxalt cos2iTnx-J2^T-C(-ak). (50) 



n=l" u k=l 
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The above expression can be checked immediately. In fact, for a = 2, i.e. (3 = 1, one has 
((—2k) = and the additional term is nonvanishing, i.e. 

A 2 (t) = 2 V / dare - * *cos27ma; = J- Ve"~. (51) 

n=l J ' V * n=l 

As a consequence, eq. fl4~6|) becomes a well known Jacobi's theta function identity. 

Furthermore, if a < 2, it turns out that A a (t) can be made very small, and we can 
write 

2£ / rfxe- |xQ|i cos27rra ~ ^ ^— ^-C(-afc), (52) 

which actually looks quite non trivial. This sum must be cut after its minimal term. 
Again, for a = 1, one can perform exactly the elementary integrations involved, and the 
result is 

£irpw=-£(*=i)i<<i-*>- (53) 

The sum on the left hand-side can be done with elementary methods and we end up with 

1 1 00 ±2r 

-cot(|) - - = _m-2r) (54) 

which is a well known but certainly non trivial identity. 

All this goes through for a = 2, 4, 6, . . .. A connection with the theory of Brownian 
processes may be stablished at that point ||. The instabilities that are known to appear 
outside the range < a < 1 and outside these particular values of a can be traced back 
in our procedure to the fact that both the full series in £'s and the additional term on the 
contour at infinity diverge. Then, both are useless, in practice, and we must come back 



to our new method as stated before (end of sect. 2) and reflected in (0). This will be 
illustrated in the following sections in great detail. 



6 Numerical estimates of the additional term 

In general, the additional term, that is, the contribution of the semicircumference 'at 
infinity' — whose existence was discovered by one of the authors a few years ago |J — has 
been very difficult to handle, even when computed at 'finite distance' Re z = c or for 
finite radius \z\ = R, either analytically or numerically (see the considerations by Actor 
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in ref. ||). Being such an elusive term to any kind of treatment, let us (for the moment) 
consider it in its most simple form 

A a = -L f dzT(z)((az) } (55) 
2m jk 

where K is the semicircumference of (finite) radius R, z = Re ie , with 9 going from 9 = tt/2 
to 9 = 37r/2. To begin with, this expression is not defined at 9 = n and one must start 
the calculation by using the well known reflection formulas for T and (. That yields 

• f dz sin(7raz/2) IYI - az) . . , . . 

A a = i — — - v . ' 1 — — — -^((l-az). (56) 

JKz(2-K) 1 - az sin(vrz) T(-z) sv ' y ' 

Such quantity is, in principle, well defined but, on the other hand, rather bad behaved. 
In modulus, the integrand clearly diverges for \z\ oo but, due to the highly oscillating 
sinus factors the final result turns out to be finite. Actually, it has been proven rigorously 
in p|, that it has the very specific value 



oo 



A 2 = -v^5(tt 2 ), S(t) = J2e^\ (57) 

n=l 

for a = 2, and that it can be made < 1 for any value of a > 2 and reasonably high values 
of R (the ones actually needed for practical applications). A possible way to handle 
expression (EST) is to employ the integral equation 



r(z)CWz) = / Sa{t)t z - l dt, Re z > 0, S (t) = T e~ n t , (58) 

what yields 

A a = -ia / dz ) 11 \ ' / S a {{2ir) a t) t~ z ~ x dt. 59 

jk sm{irz) r(— z)T{- — z) Jo 

Using now Stirling's formula and a simple approximation for the sinus fraction we get 
immediately 



A a = I dz Va (z) I S Q ((27r/a) a t) t^^dt, (60) 

V27ra Jk 



being 

ip a {z) = exp 



{2-a)z+ [ - - - 



1 n z 4- ( <\ — 2 ) : + ( — — 1 ) 7T 1 1 m : 



+ i sgn (Imz) ^— — 1^ 7r|Re z| j . (61) 
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We see that, in fact, for a < 1, ip a {z) ~~ > 0, when \z\ — > oo, while for a = 2 it is y^- 2 ) = 1 
and A 2 = — a/ttS^^ 2 ), as anticipated before (this is nothing but the famous Jacobi theta 
function identity). On the other hand, if we substitute in (|60|) its mean value, n, for the 
function ip a , we obtain 

*. = -^S.((-)°)- (62) 



a V V a 

An analytical but approximate evaluation carried out for the particular cases a = 2N, 
with N a positive integer (these values are most simple to deal with), suggests that the 
behavior of the mean value fi in terms of a can be bounded from above by an expression 
of the kind 

e — 2 Fin 

'^^tv^ (63) 

for a large. That gives for the additional term 

which is convergent for a — ► oo. 

We have checked this analytical bound with a numerical, direct evaluation of the 
integral (|SB|) written in the following equivalent form, obtained by straightforward ma- 
nipulations, 

W2 d6 sin(7rz/2) Viz + 1) w , fl 
A a = - . . - — ; 7 ; ( ^ 7 r C(^ + 1), z = Re l \ R»l. 65 
i-Tr/2 (2tt) z+1 sin(7rz/a) T(z/a) sv ; ' ' v ; 

In particular, the case a — > oo can be easily handled 

Aoo = — / — — — — sin(7T2;/2) T(z + 1) ((z + 1), z = i?e ie , R » 1. (66) 

As discussed before, it is quite difficult to evaluate such integrals for R big. Any numerical 
procedure breaks down for R large enough. This makes the method developed in sects. 2 
and 3 of this paper even more valuable. As explained in detail there, by using this tech- 
nique we never need, in practice, to go to very large i?'s in the evaluation of the additional 
term — at the expense of just calculating a few first terms of an asymptotic expansion on 
the zeta function of the differential operator one is dealing with (the asymptotic series 
must be optimally cut in the usual way). Put in the abstract (simple) situation which 



concerns us here, we need only obtain numerical values for the expressions fl65 ) and (pq) 
when R is just of the order of 20 or 30 (what is already involved enough!), since the usual 
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asymptotic series that appear in practice yield their best result for a number of terms of 
this (or very often less) order. 

There seems to be no definite regular behavior in the numerical values of the expres- 
sions above, though all of them, in the wide range 8 < R < 30, lie in the quite narrow 
interval 0.02 < |A| < 0.07. This is in very good agreement with the results obtained from 



the analytical approximation fl64|), which are A 4 ~ 0.04, A 6 ~ 0.07 and Aoo ~ 0.13. Only 
the tendency of this formula for large a deviates from the numerical result by a factor 
of 10 — namely the numerical results are about ten times lower than the upper bound 
values given by the formula QB4]). But this is not strange, since for big a's this expression 
( |64"D is rather bad (S a is then a very slowly convergent series), while, on the contrary, the 
limiting expression for A^ (|66|) remains in very good shape. 

Summing up, eqs. (|65D , (|66p and fl64|) are the best we could obtain, from an analytical 
point of view, for the treatment of the additional term which appears both in the series 
commutation procedure and in the Mellin transform method. However, for the numerical 
treatment, a clever use of the methods of sect. 3 and 4 drastically reduces the problem, 
as we have just seen in the present section: the calculation must be performed not at 



R = oo but at reasonable values of R, and in this case (p5|) and (gg) can be handled 
through standard numerical integration procedures. 



7 Example: exact summation of the string partition 
function for different ranges of temperature 

The problem of deciding if QCD is actually (some limit of) a theory of extended objects 
(strings, membranes or p-branes in general) is almost twenty years old and goes back to 
't Hooft. In a very recent contribution, Polchinski has tried to answer this question both 
for the Nambu-Goto and for the rigid strings by calculating their partition functions and 
seeing if they actually match with the one corresponding to QCD for different ranges of 
temperature [|K|. This is no place to go into the details of the procedure [fToQ , which 



will be given in a forthcoming paper and we shall here concentrate in the mathematical 
aspects of the problem only. 
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The first two terms in the loop expansion 



S, 



eff 



S + S X + 



(67) 



of the effective action corresponding to the rigid string 
1 



S 



2«n 



d 2 a 



p-i&XWXp + \ ab (d a X^d b X^ - P 8 ab ) + p / d 2 a p, (68 



being ao the dimensionless, asymptotically free coupling, p the intrinsic metric, po the 
explicit string tension (important at low energy) and X ab the Lagrange multipliers, are 
given, in the world sheet < a 1 < L and < a 2 < (3t (an annulus of modulus t), by 

L/3t 



Sn 



2a c 



A ll + X 22 f -2 + p ( 2a()/io _ X a 



(69) 



(tree level) and 



Si 



d-2 
2 

d-2 



lndet (d 4 - p\ ab d a d b 
+°° dk 



L E 



2tt 



In 



k 2 + 



An 2 n 2 
(3H 2 



4ir 2 n 2 
(3H 2 \ 22 



+ p[X n k 2 + 



(70) 



(one loop). This is a highly non-trivial calculation, which has been performed in ref. JlO 



only in the very strict limits T — > and T —>■ oo around some extremizing configuration, 
the parameters being p, A 11 , A 22 and £. As a first step of the zeta function method we 
write 



Si = -(d-2) 



where 



CaW2)=£ E [ + ™^(k 2 + y 2 + )-°/ 2 (k 2 + y 2 _r " 

_J-oo Z7T 



a 



2 \ 11 



2 P* A 

n 4- - 

+ 2a 2 



11 



A 22 V ± 



pt 



! A ll2 \ 1/2 ' 



4a 2 



1/2 



(71) 



One may consider two different approximations of overlapping validity flO|| , one for low 
temperature, (3~ 2 « p , and the other for high temperature, (3~ 2 » a p . Both these 
approximations (overlapping included) can be obtained from the exact expression above. 
It can be written in the form 



L T(s-l/2) 



2^ r( s ) 



E 



y- 



(y+y- 



F(s/2, 1/2; s; 



y- 

v = —- 

y+ 



(73) 
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This is an exact formula. For high temperature, the ordinary expansion of the confluent 
hypergeometric function F is in order 



F(s/2, 1/2; s; 1 - V ) = £ (1 - V ) k 



k=0 



k\(s) 



l + r] 



-1/2 



0. 



(74) 



since y± can be written as 
a 



V± 



t 



(n ± by + c 2 



1/2 



0t 



(75) 



After the appropriate analytic continuation, the derivative of the zeta function yields 



d_ 

ds 



Ca(«/2) 



L 



2txL 



-4^^ + fit 



1 

6 



ln(-6 2 ) + ^(-l/2)+ 7 c 



(76) 



In order to obtain this result, which comes from elementary Hurwitz zeta functions, 
C(±l, b), we have used the binomial expansion in ( [75]) what is completely consistent with 
the approximation (notice the extra terms coming from the contribution of the pole of 
to the derivative of ( A (s/2)). 
The low temperature case is more involved and now the methods of sect. 3 prove to 
be very useful. The term n = must be treated separately. It gives 



c 



(i=0), 



V2) = 77- 



L r((l - s)/2)T(s - 1/2) 



11 ,1/2-s 



2tt r(s/2) 
This is again an exact result, which yields 



s=0 



(A'V) 



2 v p 



-<Re( S )<l. (77) 



(7* 



and 



Sr° = (d-2)-y/\p. 



(79) 



Such contribution must be added to the one coming from the remaining terms (eq. 
above) 

L r(s-l/2) 



Ca(s/2) 



E 



V- 



F(s/2, 1/2; s; 1 — 77), Re(s)>l. 



0) 



v 7 ^ r(s) ^ 

(The prime is no derivative, it just means that the term n = is absent here.) Within 
this approximation, and working around the classical, T = solution: A 11 = A 22 = a fj, , 
p = t~ 2 = 1, we obtain 



C A (s/2) = -= 



L r(s-l/2) fay 



A r(s) \t 



-2s 



[F {s) + F 1 {s) + F 2 (s)] 



^1) 
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where 



Fo(s) = E 



n 



l-s 



n=l 

oo 



(n 2 + a 2 ) s l 2 



F(s/2,l/2;s;a 2 /(n 2 + a 2 ))-l 



a 



An 2 



\pt 2 p 2 
4tt 2 



a' 



n 



l-s 



Fo(s) = 



a 



2 oo 



E 



n 



-l- 



2) 



^ {n 2 + a 2 )*/ 2 ' W 4 ^ (n 2 + a 2 ) 5 / 2 ' 

The study of these functions is done in appendix A. It involves the procedures of sect. 
3 and is very related with example 2 there. Substituting the results of the appendix back 
into (PH), we obtain 



Ca(*/2) = -s 



L r(s-l/2) /a 



l-2s 



a 2 a 2 a 2 



r( s ) w 

I 1 /-co . 

— + — / dr In ( 1 - e 2 " r 
24 2?r A V 



8s 8 4 bl 

r 



— (T 



2 

+ 0(s 



a 



therefore 



Ca(0) 



Laa 2 
4t 



and 



s=0 



+ 7 



La 

T 

a 2 



a 2 a 2 a 2 



V(-1/2) T - ln(a/t)y - T - y ln(a/2) 



a 



4 12 2vr 



J dr In (l - e 27rr 



yr 2 — a 2 



Finally 



(n#0) 



(d-2 



La 

T 



lK-l/2) Y - ln(a/t)y - ^ - y ln(a/2) - | 



1 1 
+ 7 1 , 

' 4 27T J a 



dr In ( 1 



8) 



54) 



55) 



i/r 2 — cr 2 

Adding (ff9|) and (|86f) we obtain the desired expansion of the one loop effective action, Si 
near T = 0. The physical implications of these results will be discussed elsewhere. 



6) 



8 Discussion and conclusions 

The series commutation techniques that are essential in the proof of the zeta function 
regularization theorem — which, on its turn, is the basic tool in the general procedure of 
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zeta function regularization when the spectrum of the operator is explicitly known — have 
been promoted in this paper to an extremely elegant mathematical method, by making 
use of the Mellin transforms of convenient heat kernel operators, in combination with a 
rigorous treatment of the asymptotic series involved. The laborious analysis of the series 
to be commuted, the artificial picking up of a convenient function in order to mimick 
such series through pole residues on the complex plane, and the process of commutation 
itself, with the appearance of additional terms 'at infinity', have now (almost completly) 
disappeared, in favor of a quite natural Mellin transform (better, inverse Mellin transform) 
analysis of the heat kernels (see also M). Moreover, the identification of the three different 
contributions to the final result ||, namely, the naive commuted series (which results in 
a sum of Riemann or Hurwitz zeta functions), the ordinary commutation remnants (a 
polynomial function) and the very elusive, additional term of negative power-like behavior, 
appears now in a clear, natural, almost misterious way. The last contribution ||, which 
was originally very difficult to handle from a numerical point of view, has been given 
here a completely new treatment, which allows the calculation of explicit numbers with 
reasonable ease. 

Our new method has the additional advantage that it is equally well fitted for the treat- 
ment of general elliptic differential operators whose spectrum is not known (notice that 
the procedures in |2], [3], [J depend heavily on the explicit knowledge of the full spectrum). 
The sum over eigenvalues is then naturally replaced — within the same procedure — by a 
sum over heat-kernel or Seeley-De Witt coefficients. A huge mathematical industry has 
been generated for the calculation of these coefficients, and one can now get full profit 
from these result in the new context of the zeta function regularization theorem. 

The physical applications one can envisage for the techniques here developed keep 
growing every day. Aside from the many calculations, carried out in different contexts, of 
the vacuum energy and the Casimir effect in quantum field theory, condensed matter and 
solid state physics |Tl|], in the last section we have already hinted at the straightforward 
use one can make of our expressions in the calculation of the partition functions of strings 
and membranes. As pointed out by Polchinski et al. |L0 |, these results are essential in 
order to decide, once and forever, if QCD can possibly be a certain limit of some string 
(or membrane, or p-brane) theory, through the analysis of the high temperature behavior 
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of the corresponding partition functions. These kind of calculations are considered to be 
very difficult even by the mathematical-physics community (see the specific description 
in |10l|). As we have proven in the last section of this paper, they can be handled in a 
relatively simple way through our methods. We hope to be able to elaborate further on 
these and other applications in the near future. 
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A Appendix: Auxiliar expressions for the string par- 
tition function at low temperature 



We make here a carefull study of the functions Fi defined in (j82"|). F\ and F2 can be 
analytically continued without problems: 



err 

F 1 (,)=C(2 S -1)-— C(2 + l) + 5> x - 

A n=l 



y + a 2 ) -s/2_ n -s + ^_ n -s-2 



so that, in particular, 



and 



F 2 (s) = °- ((2 + 1) + £ n" 1 - (n 2 + - n 



a 



2 00 



4 



,2 J_ „2\-s/2 



which has a pole at s — 0, 



Concerning F , by using 



n=l 



limF( S /2,l/2; S ;»7(^ + ^)) = - 



l + \ 1 + 



a 



n 



17) 



19) 



(90) 



(91) 
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we can write 



1 oo 
Z n=l L 



This yields 



n=l 



- lim 

n=l . 



(n 2 + <7 2 ) 1/2 



n 



a 

2^2 



12 



+ 0(1), 



(92) 



(93) 



for the middle term. For the other two we shall make explicit use of the results obtained 
in sect. 3 of this paper. For the last term, we get 

OO 



Ef 



n=l 



n 



Res z=0 T- Z r(z)((z + 1) +C(1) = -lnr + C(l) 



(94) 



In order to apply this procedure to the first term we need some more specific knowledge 
of the function 



C-i/2W = En"*(n 2 + ^ 1/2 



(95) 



n=l 



which appears when considering (see sect. 3) 

OO OO 1 

E(^ 2 + ^) 1/2 /(™) = ^(n' + ^V'-L I dz(nr)- z M[f](z) 

n=l n=l 1711 



±-JdzT- z M[f]{z)C_ l/2 {z). 



(96) 



The study of this function can be reduced to example 2 of sect. 3. For Re z > —2p we get 
an analytical continuation which is a meromorphic function with poles at z — 2(1 — k) of 
residues ( x { 2 ) a 2k , k = 0, 1, 2, 3, . . . 



2A; 



(7 



2 A- 



fe0V ,y, + 2 ( ,-T) + ^ 2 + CT2 » 1/2 -r( £2 + CT2 » 1/2 

i ?L ) *ln(l-e—) d 



27T *j?JC<.±) 



dt 



(97) 



where the three last terms are integer functions and the contours are chosen avoiding 
the points ±ia. 

In our case it is sufficient to take p — 1. Considering now — 2 < Re z < 0, the limit 
e — > can be taken naively, with the result 



C- 1/2 (z) = j o drr 1 '* 



i + ^-f(p-i)E 



fe=0 



l/2\ fa 
k , 



21: 



(98) 
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E 



1/2* 



(7 



2A- 



— - — V / dt In fl - e ±2mt ) 4- \t~ z (t 2 + a 2 ) 1 ' 2 



fe 4V fc / * + 2(*-i) ± 

It is the last term the one which prevents the continuation to the right of z = due to 
its divergence at t = 0. This needs a special treatment (see again sect. 3), and the final 
result is 



„27rr 



\Jr 2 — 



+ (99) 



The integral term is awkward but harmless: it is exponentially suppressed as compared 



to the rest. Now we can go back to (|9q) 



E(- 2 + 



a 



2\l/2 p -rn 



1/2 W 



n=l 



7<x 2 cr 2 a 2 , /cr 

■- 1 In - 

2 4 2 V2 



1 z" 00 / 
H — / dr In ( 1 — e 

7T Jcr V 



1 °\ 

--y lnr 



a 
2 



r 

2nr 



(100) 



r 



\/r 2 — cr 2 



+ a(i). 



And this yields for F (0): 

r 2 



^0(0) 



7<T 



a 2 a 2 , 
H In 

4 8 4 



a 1 1 r°° 
4 24 + 2 



- / dr In ( 1 - e 2 " r 

7T Jcr V 



\/r 2 — a 2 



(101) 
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